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ABSTRACT 

Based  on  the  wavelet  transformation  and  adaptive  Wiener  Filtering,  a  new  method  was  presented  by 
the  authors  to  perform  the  synchronization  and  detection  of  the  binary  data  from  the  Free-Space 
Optical  (FSO)  signal  [1],  It  was  shown  in  [1]  that  the  Haar  wavelet  with  a  fixed  scale  is  an 
excellent  choice  for  this  purpose.  The  output  of  the  filter  is  zero-mean  and  is  closely  related  to  the 
derivative  of  the  binary  data.  In  this  effort  an  analysis  of  the  work  in  [1]  is  presented  to  obtain  the 
probability  of  bit  error  using  a  Bayesian  ternary  hypotheses  testing.  The  analysis  also  results  in 
determining  optimum  thresholds  for  the  detection  of  binary  data. 
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1.  INTRODUCTION 

Free-Space  Optical  (FSO)  communication  system  is  an  important  area  of  research  due  to  its 
important  advantages  of  providing  a  very  large  bandwidth  and  relatively  low  cost  of 
implementation.  These  systems  have  many  desirable  applications,  perhaps  the  most  important  of 
which  is  their  use  in  providing  connection  from  the  high  bandwidth  fiber-optic  backbone  to  the 
buildings  and  businesses  desiring  high  bandwidth  access  without  the  high  cost  of  installing  fiber 
through  the  local  infrastructure.  The  FSO  communication  systems  also  have  numerous  advantages 
over  the  Free-Space  RF  (FSRF)  systems. 

Bit  rates  of  FSO  systems  are  approaching  100  Gbits/sec  [2-5].  Besides,  atmospheric  propagation 
effects  cause  the  FSO  channel  to  be  ‘bursty’  in  nature  and  highly  variable  with  a  rate  of  change  as 
high  as  1  kHz  and  with  power  fades  greater  than  10  dB.  Therefore,  because  of  a  very  high  bit  rate 
and  a  time-varying  channel,  the  bit  error  reduction  is  a  difficult  and  vital  area  of  research.  The 
purpose  of  this  paper  is  to  introduce  and  implement  a  new  and  effective  method  to  deal  with  these 
issues. 
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In  a  typical  RF  digital  communication  problem,  the  binary  data  is  designed  to  avoid  Inter-symbol 
Interference  (ISI)  and  is  digitally  modulated.  Here,  the  channel  is,  in  general,  modeled  as  band- 
limited,  Additive,  White  and  Gaussian  Noise  (AWGN),  and  the  demodulated  data  consists  of  a 
string  of  pulses  with  a  constant  power.  The  value  of  the  mean  is  typically  zero  and  the  variance  of 
the  data  is  a  constant.  The  detection  problem,  therefore,  is  to  use  matched  filtering  and  constant 
thresholding  to  identify  the  digital  data.  In  contrast,  the  FSO  communication  channels  can  be 
modeled  as  wide-band  and  time-varying  systems.  The  selected  format  of  the  binary  data  is  typically 
NRZ-L“  [6,  p294],  with  no  sinusoidal  modulation.  The  received  FSO  signal  consists  of  a  high 
frequency  data,  which  is  varying  in  its  amplitude.  This  corresponds  to  a  non-stationary  data  with 
variable  and  irregular  values  of  mean  and  variance.  Consequently,  the  conventional  matched 
filtering  and  constant  thresholding  is  not  applicable.  Figure  1  displays  a  typical  received  FSO 
signal.  This  figure  shows  that  the  binary  data  is  a  noise-like  signal  with  variable  variance  and  mean. 
The  data  can  be  viewed  as  a  high  frequency  signal  modulating  a  relatively  much  lower  frequency 
waveform.  In  this  paper  an  analysis  of  the  work  in  [1]  is  presented  to  obtain  the  probability  of  bit 
error  using  a  Bayesian  ternary  hypotheses  testing.  The  analysis  also  results  in  determining  optimum 
thresholds  for  the  detection  of  binary  data. 


2.  MATHEMATICAL  MODELING  AND  PROCESSING  OF  FSO  SIGNAL 

The  received  FSO  data  can  be  mathematically  described  as  follows: 

r{t )  =  <7(t)d  (, t )  +  mL  (t)  +  mH  (t)  (1) 


where  d(t )  represents  the  binary  data  (message)  and  a (t)  signifies  the  variation  of  the  signal 
amplitude  due  to  atmospheric  transmission  effects.  In  addition,  the  model  assumes  two  types  of 
additive  noise.  The  first  noise,  mL(t) ,  stands  for  the  relatively  low-frequency  fluctuations  of  the 


signal  mean  value  caused  by  insufficient  ac  coupling.  The  second  term,  mH(t) ,  characterizes  the 
additive  white  Gaussian  noise  (AWGN)  with  zero-mean.  As  reported  by  the  authors  in  [1],  and 
depicted  in  Figure  2,  we  first  pass  the  data  r(t)  through  an  adaptive  Wiener  Filter  (AWF)  for 
denoising.  The  output  of  this  filter  is  then  passed  through  a  Haar  Wavelet  transformer  with  the 
transfer  function  [1]: 


gjaO)!  2 


(2) 


‘One”  is  represented  by  level  “+1”  and  “Zero”  is  represented  by  level  “-F 
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where  a  is  the  scale  of  the  Haar  wavelet.  Figure  3  depicts  a  typical  FSO  data  at  the  input  and  the 
output  of  the  filter  (2)  with  fixed  scale  a.  As  shown  in  this  figure,  the  filter,  shown  by  (2), 
efficiently  removes  mL(t)  in  (1),  and  transforms  the  received  data  into  a  zero-mean  signal 
resembling  an  approximation  of  its  derivative.  The  filtered  data,  as  depicted  in  Figure  3,  is  then 
sampled  at  the  multiple  integers  of  bit  time,  assuming  that  it  has  been  properly  synchronized  [1],  In 
the  following  section,  we  establish  a  ternary  hypotheses  test  to  calculate  the  probability  of  bit  error. 


3.  CALCULATION  OF  PROBABILITY  OF  BIT  ERROR 

This  section  is  devoted  to  the  calculation  of  the  probability  of  bit  error.  In  view  of  Figure  3,  the 
output  y  of  the  sampler  shown  in  Figure  2  can  be  considered  as  a  scalar  random  variable,  which  can 
be  characterized  by  a  ternary  hypotheses  test  as  follows3: 

H0:  y  =  w, 

Hx\  y  =  +d.f1+w,  (3) 

H2:  y  =  —d.f2  +w, 

where  (3)  is  followed  from  (1)  assuming  that  the  residual  noise  at  the  output  of  the  sampler  is 
characterized  by  w  and  is  considered  to  be  white  and  Gaussian  with  mean  zero.  In  addition,  f] , 
i  =  1,2  are  Rayleigh  fading  random  variables,  independent  of  w,  and  d  represents  a  known,  non¬ 
negative  quantity  proportional  to  energy  of  the  binary  data.  Without  loss  of  generality,  we  assume 
that  d  —  1 .  Hypothesis  H0  corresponds  to  a  scenario  for  which,  at  the  transition  time,  there  has  not 

been  any  bit  change  (d  =0) ,  and  hypotheses  H]  and  H2 ,  on  the  other  hand,  represent 
circumstances  for  which  bit  transitions  have  occurred  from  low-to-high  and  high-to-low, 
respectively.  Based  on  the  aforementioned  assumptions,  we  can  establish  a  ternary  hypotheses  test 
by  determining  the  conditional  probability  density  functions  (pdf)  of  y  given  the  three  hypotheses 
H0 ,  H]  and  H2  [7  P23].  A  maximum  likelihood  detector  can  then  be  implemented  by  evaluating 
these  pdfs  at  any  given  observation  Y  and  select  the  largest. 


3  We  use  lower  case  letters  for  a  random  variable.  Upper  case  letters  are  utilized  to  signify  the  running  values  of  the 
corresponding  random  variable.  Also,  the  model  (3)  considers  the  fact  that  the  fading  channel  may  act  differently  for 
high  and  low  bits. 
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To  begin,  note  that 


Accordingly, 


Py(Y\H0 ) 


Py{Y\HvF,) 


1  „-(Y-Fi)2 12<rl 


V2^2 


py(Y\H2,F2)  = 


-(Y+F2)2/2crl 


Also,  by  definition,  /'. ,  i  =  1, 2  are  Rayleigh  distributed  random  variables;  therefore, 


pf(Fi)=-^ie  u^Fi) 


where  u{.)  is  the  step-function.  It  follows  from  (5)  and  (7)  that 


(Y-F^fnal 


pf(F1)dFl 


where 


_2  .  2 

_cjh.  +  ^/1 

(T»(T/i 


<2<A)  =  1 —  fe  1  l2dx 

.  /O  -77"  J 


(10-a) 


with  the  known  property  that 
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Q(x)  +  Q(-x)  =  l 


(10-b) 


Similarly,  it  can  be  shown  that 


**  'H^  p’iY  1  M 


wherein 


_2  .  _2 
<7  +  <7  f 

n  w  J  2 

^2  ~  2  2 

a«°h 


It  is  noted  that  the  value  of  Y  for  which  py(Y  \  Hx)-  py(Y  I  H0)  characterizes  the  threshold  value 
7j .  Similarly,  the  values  of  Y  for  which  p  (Y  I  H2)  =  p  (Y I  H0)  signifies  the  threshold  value  T2. 
Consequently,  it  follows  from  (8)  that  at  Y  =  7j  we  attain 

(i3) 

[  4Bx(7~w  jBt(7w  J 


Correspondingly,  it  follows  from  (11)  that  at  Y  =  T2  we  obtain 


Unu 


l  \  Jbi<£, 


)  =i 


For  convenience,  define  an  auxiliary  function  gj(Z)  as  follows: 


g1(Z)  =  Al{l  +  V2^Zcz2/22(-Z)} 


1  Y 

where  A  = - —  and  Z  =  — == — .  It  follows  from  (15)  that  the  value  of  Z  =  £  for  which 

g,  (Z)  =  1  can  be  used  to  determine  the  threshold  value  7j ;  that  is, 


Note  also  from  (9)  and  (16)  that 
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in  which  y/  represents  the  value  of  Z  for  which  g2(Z )  =  1 . 

Once  Tt  and  T2  are  calculated,  we  can  obtain  the  average  probability  of  error  which  can  be 
obtained  by  considering  the  conditional  probability  of  incorrect  reception  P(e  I  Ht)  as 

P(s)  =  ^[P(e  I  H0)  +  P(s  \Hl)  +  P(£\H2)]  (20) 

where,  without  loss  of  generality,  it  is  assumed  that  the  hypotheses  are  all  equally  likely.  It  follows 
from  Figure  (4),  and  (10-b),  that 
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observation  Y 


Figure  4.  Conditional  pdf’s  of  the  Ternary  Hypotheses 


P(£\  H0)  =  1— 7=!=^  ['e^'^dY  =  l-Q(T2 /(Jj  +  QiT,  /  aw) 
=  Q(-T2/crJ  +  Q(Tl/(jw) 

Also,  in  view  of  Figure  4,  it  follows  that 


P(e\Hl)=^py(Y\Hl)dY 


1 


Py(Y\H0 ) 


/2^<rJ 


In  Appendix  A,  it  is  shown  that  (22)  becomes 


(21) 


(22) 
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P(s\H1)  = 


_2  .  _2 
cr«'+cr/1 


_2  .  _2 

CrH’+Cr/1 


o--;  ~  A^c\  n 

-^-T-  V  1  /lJ  1- 2(f)  - 
^  +  (Ta  a*  +  ah 


r  (o'1 + oi )  i 

1-G  ,  ;i  <T 


Using  (10-b),  (23)  can  be  simplified  as  follows: 


1  l  +  (o2f/o2w)) 

P(elH,)=  — — — -  \q  -  V'  ;  ( 

I  (V'U  J 


(°F<) 

1+K/cr») 


_I  £«L  U 

e  2UJ  G(_^)_ 


1+K/o0 


1+K/g0 


It  is  noted  from  Figure  4  that 


P(e\H2)=\py(Y\H2)dY  =  \Py(-Y\H2)dY. 


It  is  eminent  from  (8)  and  (11)  that  py(-Y\H2)  and  py(Y  I  Ht )  have  the  same  forms  when  we 
replace  B]  with  B2 ,  and  <72  with  <72.  .  Hence,  in  view  of  (17),  (19),  and  (25),  P(s\H2 )  can  be 
realized  from  P(e  I  Hx) ,  shown  by  (24),  by  switching  Bx  to  Bn ,  o2f  to  <j2f  ,  and  g  to  -y/ ;  that  is, 


1  ]  f  1  +  (cr2  /cr) 

P(£\H2)=  — — - —  \q  /  f — p- 

1+(<t/2/<t„)|  I  (ollcl) 


V  - 


K/g.) 

1+K/<jt) 


-MV  I"  {o2lo2\  V 

\e  ^  Q(Y)~  }/\  *L  ! 


1  +  K7^) 


1  +  K/c7^) 
K7^) 


Hence,  the  total  average  probability  of  error  can  be  calculated  by  substituting  (21),  (24)  and  (26)  in 
(20).  Finally,  it  should  be  mentioned  that  for  the  case  for  which  <72.  =  a)  =  o2f ,  it  follows  that  Bx 
=  B 2  and  7j  =  —T2  =  T  .  Consequently,  it  is  seen  from  (20)  and  (21)  that  P(e)  becomes 
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(27) 


P(£)  = 


Q(—)  +  P(£\H1 ) 
cr 


In  view  of  hypotheses  Hx  and  H2  in  (3),  we  define  the  Distorted  Signal-to-Noise  Ratio  (DSNR)  as: 


_2  2 

(7,(7 , 

DSNR.  =  -^—A- ;  i  =  1,2. 


(28) 


Since,  by  definition,  d  =  1 ,  it  follows  from  hypotheses  //  and  H2  in  (3)  that  crj  =1.  Consequently, 


from  (28)  we  have 


DSNR.t  = 


i  =  1,2. 


(29) 


For  convenience,  we  assume  that  <J2f.  =  <7 j  =  cr2f ;  hence,  DSNR]  =  DSNR2  =  DSNR.  Figure  5 
illustrates  the  total  average  probability  of  error  (27)  as  a  function  of  DSNR .  This  figure  shows  that 
as  <72  /  cr  increases,  the  probability  of  error  decreases.  Figure  6  depicts  the  normalized  threshold 

T / o  as  a  function  of  <7]  la1.  It  is  seen  from  this  figure  that  the  normalized  threshold  T / <7  is  a 

XV  J  w  c?  w 

2  2 

monotonic  function  of  cry  /  cry .  A  practical  implication  of  this  result  is  that  the  optimum  threshold 
can  be  evaluated  from  Figure  6  by  a  table-look-up  strategy.  That  is,  once  cr 2  /  <7 1  is  estimated,  the 
normalized  threshold  is  found  from  Figure  6,  and  the  optimum  threshold  is  estimated  by  multiplying 
the  result  by  crw. 
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Figure  5.  Average  Probability  of  Bit  Error  Versus  DSNR 
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DSNR 


Figure  6.  Normalized  Optimum  Threshold  Versus  DSNR 
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Appendix  A  -  Derivation  of  P(e  I  Ht ) 


The  objective  of  this  appendix  is  to  illustrate  the  derivation  of  P{e  I  H t)  shown  in  (23).  We  know 
from  (8)  that 


1 


py(Y  I  Ht)  =  — -py(Y  I  H0)V  +  -^Yerl/2*a4"Q(-r?-T) 

5,(7“ 


V2 ~K 


\[Bxc 


(A-l) 


Hence,  the  conditional  probability  of  error  P(e  I  Hx )  becomes 


1  71 

P(£\Hl)  =  P  py(Y  \  H^Y  =  — — —  f  PV(Y  I  H0)dY 

By°h  i 

. 


(A-2) 


Application  of  (4)  in  (A-2)  renders, 


P(£\H.)  =  -^[l-Q(-^-)]  +  /c  f  e~a  Y  YdY  -  K  [e^YQ^YMY 

<TW  i  i 


wherein 


(A-3) 


a 2  = 


K 


1  1 


2cr  25,  a 

w  1  w 

1 


P 


1 


The  last  integral  in  (A-3)  can  be  solved  using  integration  by  parts;  that  is, 

\udv  =  uv-  j  vdu 


with 


m  =  Q(fiY) 
dv  =  e~a  Y  YdY. 


(A-4) 


(A-5) 


(A-6) 


Hence, 
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}  e-a2rlYQ(j8Y)dY  =  -^Q^T,)e^  -  )^-2 


**  f  1  e^e^dY. 


(A-7) 


Furthermore  the  second  integral  in  (A-3)  is  solved  as 


K  -a2T2 


K  e~a  YdY  = - -e 

J  2a 


(A-8) 


Applications  of  (17),  (A-4),  (A-7),  and  (A-8)  in  (A-3),  after  collection  of  terms  and  further 
simplification,  finally  leads  to  (23). 
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